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Abstract. Second-order dissipative hydrodynamic equations for each component of a multi-component 
system are derived using the entropy principle. Comparison of the solutions with kinetic transport results 
demonstrates validity of the obtained equations. We demonstrate how the shear viscosity of the total 
system can be calculated in terms of the involved cross sections and partial densities. Presence of the inter- 
species interactions leads to a characteristic time-dependence of the shear viscosity of the mixture, which 
also means that the shear viscosity of a mixture cannot be calculated using the Green-Kubo formalism 
the way it has been done recently. This finding is of interest for understanding of the shear viscosity of 
a quark-gluon-plasme extracted from comparisons of hydrodynamic simulations with experimental results 
from RHIC and LHC. 
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1 Introduction 

The deconfined state of QCD matter produced at the early 
stage of ultrarelativistic heavy-ion collisions at RHIC and 
LHC is a multi-component system with quark and gluon 
degrees of freedom. Large values of elliptic flow coefficient 
«2, observed at RHIC [TJ[2] and LHC [3], indicate that 
the produced quark-gluon plasma (QGP) is a nearly per- 
fect fluid. This has motivated rapid developments on rel- 
ativistic dissipative hydrodynamic formalisms [HOEJEHS] • 
The value of the shear viscosity to the entropy density 
ratio rj/s for the QGP at RHIC and LHC was extracted 
from comparisons of hydrodynamic [§lllO| as well as ki- 
netic transport |llj calculations with experimental data. 
All these hydrodynamic formalisms are based on the as- 
sumption that the quark-gluon mixture can be regarded 
as an effective one-component system, where r]/s is an 
external parameter characterizing the dissipation in the 
system. For a one-component system the rj/s ratio can be 
calculated using the Green-Kubo formalism [THITB] as well 
as other systematic approaches |14lfT5] . One may ask the 
question whether analogous calculation of the shear vis- 
cosity of a mixture is possible, i.e. whether a mixture be- 
havior is equivalent to a one-component system or rather 
not. These questions are of major interest for investigation 
of the QGP properties. 

In this paper we demonstrate that a standard one- 
component hydrodynamic description with a single shear 
viscosity coefficient calculated by e.g. Green-Kubo formal- 
ism in general cannot be applied to a multicomponent sys- 
tem. We will explain this statement by deriving second- 
order dissipative hydrodynamic equations for a multi - 



component system from the entropy principle. Our ap- 
proach differs from the one reported in Ref. [5], since we 
introduce separate evolution equations and transport coef- 
ficients for each component of the mixture. We then show 
that by summing-up equations for all components one can 
obtain an equation for the system as a whole, which has a 
relaxation-type form characteristic for all the present hy- 
drodynamic formalisms, but the effective shear viscosity 
for the mixture is now related to the partial shear pres- 
sures of its components and thus has a non-trivial time 
dependence which is not supported by the Green-Kubo 
or, for this matter, any other formalism, in which an equi- 
librium state of matter is assumed. We will confirm our 
findings by calculating the shear viscosity of a mixture us- 
ing the Green-Kubo formula and comparing solutions of 
the hydrodynamic equations with the ones from kinetic 
transport calculations. 



2 Dissipative hydrodynamic formalism for a 
multi -component system 

We consider a mixture of N particle species, for which we 
define a common velocity field u^ L . Neglecting bulk pres- 
sure and heat flow we can construct the total entropy cur- 
rent as mi 
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e|t^<v, (i) 



where s eq is the total entropy density in local equilib- 
rium and is the hydrodynamic velocity. Tj and are 
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the temperature and local energy density of the particle 
species i. In analogy to the one-component case [51IT5] one 
obtains ft — (9/4ej). tt^ = T? v — T^" q is the shear stress 
tensor, which is, as long as heat flow and bulk pressure are 
neglected, the difference between the energy-momentum 
tensor T^ 1 and the equilibrium one. Equation ([1]) is the 
generalization of the entropy current for a one-component 
system (TV = 1), discussed for instance in Refs. [^1151115] . 
The total entropy production is then 



a<*f> 



V 2T 4 J Ti M 1 



with the shear tensor 



(2) 
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and A a /3 = g a f3~u a up with the metric g a p — diag(l, —1,-1, 
We have used conservation of the partial particle flows and 
total energy-momentum tensor, d^Nlf — and d^T^" = 
0, to obtain Eq. ©. The obtained equation for the entropy 
balance is again a generalization of the one-component re- 
sult 0. 

According to the second law of thermodynamics the 
entropy production is non-negative. A simple way to fulfill 
this is to make the terms in the square bracket in Eq. ((2j 
to be proportional to 7r"^, [•••] = nf 13 j (2r)iTi), Then the 
entropy production © has the following algebraic struc- 
ture: 

>0. (3) 
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This leads to the dynamical equation for each 7r, 



a/3. 
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which is analogous to the equation introduced by Israel 
and Stewart for a one-component system (N = 1) [HE]- 

In order to apply Eq. ||J to a multi-component sys- 
tem (N > 1), we first need to determine the coefficients 
rji, which in general differ from the usual definition of the 
shear viscosity. The reason for this is that 7r"^'s are cor- 
related due to interactions between particles from differ- 
ent species. These correlations between irf s can only be 
seen, when each rji depends on all tt" , j = 1, 2, • • • , N. 
We will also show later that the coefficients rji become the 
shear viscosities, when the ratios of components of n"^'s 
are relaxing to constants in time. 

We now make use of relativistic kinetic theory and ex- 
press the entropy current via the phase-space distribution 
function fi(x,pi): 

N 

s " = E / dr * rfM*>Pi) I 1 - tefi&Pi)} (5) 

i—l 

with ATi = d 3 pi/ Ei/ (2ir) 3 . It was shown for the case of 
N = 1 [T7] and is obviously true for N > 1 that using the 



Grad's ansatz [5] 

fi(x,Pi) = f i<eq (x,pi)(l + AiKi^upfp") (6) 

in Eq. ([5]) one obtains Eq. (TTJ) up to second order in i^i,^. 
Here ji fiq (x,pi) is the equilibrium distribution function 
and Ai — [2(ej -I- P^Tf] -1 [TS], where Pi is the pressure. 

The space-time evolution of fi{x,pi) obeys the Boltz- 
mann equation 

N 

tid^ = h, ■■■ j n } = Cum + J2 c * /j] > 

(7) 

where Cu are the collision terms describing interactions of 
particles of same species and describing binary interac- 
tions of particles of different species. Explicit expressions 
forxthe collision terms can be found for example in |18) . 
Taking derivative of ([5]) and using ([7]) we obtain 



N 

^ ^ A-i7Ti ^ LU 
i=l 



(8) 



Comparison between Eqs. (|S]) and (J3J leads to 



2Aim^ J dr iP t*p?Ci 



(9) 



Because the collision term C, is a functional of all /,'s, 



each r]i depends on all ttj 's with j = 1, 2, • • • , N. 

In order to simplify Eq. we will now consider a one- 
dimensional system, which implies that the shear stress 
tensor has a diagonal (and of course traceless) form: 7r Q ' 3 = 
diag(0, 7r/2, 7r/2, — 7r), with the single independent compo- 
nent 7T. By virtue of Eq. ([6]) this form of the shear stress 
tensor is equivalent to a deformation of the momentum- 
space distribution along the z axis, whereas in the trans- 
verse xy plain the momentum distribution is isotropic. 
Furthermore we consider only two species, i,j — 1,2 and 
isotropic scattering processes, i.e. we assume that the dif- 
ferential cross section dajdfl does not depend on scat- 
tering angle. Inserting the off-equilibrium destribution © 
into the collision term in (|9]) and using the aforementioned 
simplifications we obtain the following expression for the 
shear viscosities of the mixture constituents: 



N 



7 n. 
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(10) 



In the latter equation nj/ni denotes the ratio of the par- 
ticle densities of mixture constituents. In our formalism 
this chemical composition is fixed by the initial consition 
and does not change in time, although this assumption is 
a very strong simplification for a QGP. Now the obtained 
expression for r\i can be inserted into the dynamic evolu- 
tion equation Q for TTi (that is we take the ^33 component 
of 7r M „ in although considering any other component 
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will lead to identical result). We obtain a dynamic equa- 
tion for the shear tensor components 7r$ in a mixture: 



5 7 \ 2 

-niau + -nj<Tij I TTi + -ni<TijTTj 
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(11) 



with the short-hand notation tt = u^d^n and a denoting 
the zz component of the shear tensor a = er 33 . A distinct 
feature of the obtained equation, which has a characteris- 
tic relaxation-time form, is the presence of two time scales 
for relaxation, of which the second one is associated with 
a coupling between the partial shear pressures TTi of the 
two species. 



3 Coupled dynamics in a mixture 

Now the question arises, as of how the dependence of r\i 
in Eq. (|10[) on the ratio fti/TTj can be interpreted. If rji is 
understood as the shear viscosity of the medium, it must 
be a property of the medium, i.e. depend exclusively on 
its chemical composition and the associated cross sections. 
Moreover, the shear viscosity must be defined in the prox- 
imity of the equilibrium. Let us consider the following sit- 
uation: all velocity gradients vanish. This means that the 
shear tensor a^ v , and thus a as well, vanish. We also as- 
sume the the temperatures of the two sub-systems are 
equal and can be replaced by the single temperature T 
of the system, i.e. Tj = T. We move into the fluid rest 
frame, in which u^d^ = d/dr. Equations (jTTJ) are now re- 
duced to a set of two coupled differential equations for the 
shear stress tensor components m, 7r 2 with the parame- 
ters n±, n 2 , en, <7i2 and tx 22 . The system of equations can 
be solved analytically and thus a solution for ^ (r) is also 
found: 

— (t) = A(n, a) • tanh (r • B(n, a) + V(n, a, vr )) . (12) 

In the latter equation A, B and T> are algebraic func- 
tions of the system properties and o~ij , and in particular 
T> depends on the initial value of the ratio 7Ti/7r 2 . Equa- 
tion (fT2"]) demonstrates that the ratio 7Ti / 7r 2 of the partial 
shear pressures is determined completely by the properties 
of the system, but is not a constant. Moreover, Eq. (|12l) 
leads to saturation of 7ri/7r 2 in r — > 00 limit, i.e. a well- 
defined characteristic limit for the ratio exists. In general 
this characteristic value is different from the value n\ /ri2 ■ 
It is also important to mention, that there is no conclusive 
way to specify the initial values of tx\ and 7r 2 - or just 7r 
for a standard one-component hydrodynamic calculation. 
In most hydrodynamic approaches the standard choice is 
therefore the trivial initialization tt(tq) = 0. Another pos- 
sible choice is initialization with the Navier-Stokes value 
7t(to) = 2rja [TO]. As it is for the one-component case, for a 
multi-component system the choice of the initial condition 
for TTi is not clear as well. If the gradients are 'switched 
on', the trivial choice 7r,(ro) = will lead, according to 



Eq. dTTJ) to TTi = a J 'Pi. And since fc/Pi = ei/e 2 = ni/n 2 , 
after a short time the shear pressure ratio will be approx- 
imately equal to the density ratio, 7Ti/7T2 rs ni/n^- 

It is also interesting to build the sum of Eqs. (fTTj) . 
which leads us to a relaxation equation for the total shear 
pressure in the mixture t\ = tt\ + 1x2- If we write this 
equation in the relaxation-time form with only one relax- 
ation scale, as is usual for most hydrodynamic formalisms 
presently used [9lfTU]. we obtain 



5 

1 



1 + 7Tl /tt 2 



Kl/K2 
+ 7T1 j%2 

(n 2 <J22 + ni<ri2 



(ni<7ii + n 2 ai 2 ) 

- gradients . (13) 



According to the standard definition of the relaxation time 
7V in the second-oder formalisms [HI 



9r7 
2e 



(14) 



we recognize that the shear viscosity of the mixture is now 
given by 



7]mix — r 6 
5 
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l + 7ri/7r 2 l + 7Ti/7r 2 



a; 



, (15) 



with the inverse of the mean free path A^ 1 = n\0\\ +ri 2 o"i 2 
and A^ = n 2 (722 + n\0\% and the time-dependent ratio 
t^i/^2 given by Eq. (fT2")l . 

The obtained result for the shear viscosity of the mix- 
ture is interesting for the following reason. If one attempts 
to calculate the shear viscosity of the mixture we have 
considered here, e.g. using the Green-Kubo formula in a 
kinetic transport simulation, the result would naturally be 
a constant value for r) m i x . On the other hand the result 
obtained here in form of Eq. (|T5|) implies that the viscosity 
of a mixture is time-dependent, but saturates. The value 
at which the time-dependence dies off will not be identi- 
cal with the value one would obtain using the Green-Kubo 
formalism because the coupling between the species in a 
mixture induces an internal dynamics in a system. 

To demonstrate this, we calculate the shear viscosity 
coefficient of a mixture in the kinetic transport model 
Boltzmann Approach to MultiParton Scatterings (BAMPS) 
[TSlfTi"] using the procedure successfully applied by us in 
Ref . [12] • The cross sections for the scattering processes of 
the two species, confined in a static box, are chosen to be 
an = 10 GeV~ 2 , <j 12 = 5 GeV' 2 and cr 22 = 2.5 GeV- 2 . 
The density ratio is ni/n 2 = 5. The temperature is chosen 
to be T = 0.4 GeV and both particle species are consid- 
ered to be Boltzmann gases with degeneracy factors 16, 
i.e. m = 5/6 • 16/tt 2 T 3 and n 2 = 1/6 • 16/tt 2 T 3 . Equation 
of state is the ideal one, i.e. e\ — 2>ri\T and e 2 = 3n 2 T. 
Note that for this setup the mean-free path scales for the 
two species are Ai = 0.207 fm and A 2 = 0.414 fm. These 
values are chosen to crudely simulate quarks and gluons 
in a QGP. To obtain the shear viscosity, we extract the 
correlation function 

C(r) - i {{tt xv {Q)rr xy {t)) + (n xz (0)tt xz (t)) + (tt vz (0)n vz (t))) 

(16) 
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Fig. 1. Autocorrelation function extracted from BAMPS 
(symbols) as function of correlation time. Two-exponent fit 
Eq. |JIgJ with Ci = C* 2 = 8.66 ■ f0" 4 G*eV 2 / fm 6 , n = 
0.264822 frn/c and r 2 = 0.479556 fm/c is shown by bold 
solid line. The two exponents from the two-exponent fit are 
also shown separately by the thin solid and dashed lines. 



where r is the correlation time and (.) denotes ensemble- 
average in the static box. The Green-Kubo formalism re- 
lates the shear viscosity to the integral of the correlation 
function over the relaxation time: 



oo 

V=^J C{r)dr, 



(17) 



where V denotes the considered volume and T the tem- 
perature of the system. In a standard one-component case 
the correlator C(r) is very good described by an expo- 
nential function with the relaxation time r c |121I13| . For a 
mixture, however, we find that a single-exponent fit does 
not work and two-exponent fit must be considered: 



C(t) = d ■ e- T/Tl + C 2 ■ e~ T/T2 



(18) 



with two relaxation times t% and ti. In Fig. [JJ we de- 
montrate the correlation function extracted from BAMPS 
static box calculations together with the two-exponetial 
fit. Integrating the correlation function shown in Fig. [JJ 
over the correlation time we obtain for the shear viscos- 
ity i] = 0.062 GeV~ 3 . To demonstrate the difference be- 
tween this result and the effective mixture viscosity (|15p 
we show in Fig. [2] the time evolution of rj m i X for the setup 
described above. From Fig. [5] we recognize that the mix- 
ture viscosity T) m i x is approximately equal to the Green- 
Kubo result at erly evolution stage but significantly in- 
creases with time. In order to verify the applicability of 
the obtained relaxation-type equations ([TTj) for the shear 
pressure we now study relaxation of shear pressure in a 
static BAMPS box. We use a quasi-static setup, i.e. a vol- 
ume with no gradients but finite initial shear pressures 
7i"i (0) and 7^(0). In the kinetic transport solver BAMPS 
this is achieved by sampling particles isotropically in space 



Fig. 2. Time evolution of the shear viscosity r} m ix from Eq. 
(|15[) . See the text for details of the setup. Result of the Green- 
Kubo formalism is indicated by the arrow. 
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Fig. 3. Relaxation of the total and partial shear pressures in 
BAMPS (symbols) and analytic solutions (lines). Additionally, 
one-component solution with the viscosity obtained using the 
Green-Kubo formalism is shown to reproduce the total shear 
pressure of a multi-component system only at early times of 
the evolution. 



according to the distribution function (J6j) with a chosen 
value of 7r = ir zz . This setup provides a cross-check of 
the relaxation dynamics described by the first two terms 
on the right-hand-side of Eqs. (fTTj) and also provides a 
possibility to cross-check of the validity of Eq. (fT2"1) . We 
use same cross sections and composition of the gas as for 
the Green-Kubo calculations discussed above. Results of 
BAMPS static box calculations are shown in Fig. [31 The 
BAMPS results (symbols) for the total as well as the par- 
tial shear pressures (Fig. [3] (a)) are well reproduced by 
the analytic solutions of Eq. (fTTj) . The ratio tti/tt2 calcu- 
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lated in BAMPS is also shown to agree with the analytic 
solution and demonstrates the expected saturation (Fig. 
[3](b)). The one-component solution 7t(t) = 7r(r ) • e~ T / T,r 
with the shear viscosity value tjck calculated using Green- 
Kubo formalism is shown by the solid black line in Fig. [3] 
(a) and is unable to describe relaxation of the total shear 
pressure in a multi-component system on a long time scale. 
This means that the Green-Kubo formalism in the form it 
is applied to one-component systems cannot be applied to 
calculate the shear viscosity of a mixture. From Fig. [3] (a) 
we rather recognize that the viscosity must increase with 
time, though it might be close to the Green-Kubo result 
at early times of the evolution. 



4 Conclusions 

We have derived second-order hydrodynamic equations for 
the shear tensor components of constituents of a mixture. 
A cross-check of the obtained equations is provided by 
comparisons of the solutions with kinetic transport calcu- 
lations with BAMPS, which demonstrate very good agree- 
ment of the results. We have also demonstrated that the ef- 
fective shear viscosity of a mixture of two components does 
have a non-trivial time-dependence, which is explained by 
inner dynamics of the mixture due to inter-species interac- 
tions. If the Green-Kubo formalism is applied to calculate 
the shear viscosity of a mixture, the result cannot capture 
its time-dependence. Thus, if calculated by the Green- 
Kubo formalism, the shear viscsosity cannot be used to 
describe hydrodynamic evolution of a mixture. It will be 
very interesting to investigate the impact of our finding on 
extraction of the shear viscosity of a Quark-Gluon mixture 
from comparisons of experimental data with the results of 
dissipative hydrodynamic simulations. We expect for ex- 
ample, that the elliptic flow coefficient V2(pt) calculated 
from a mixture, such as the Quark-Gluon plasma, can- 
not be accurately reproduced by a one-component dissi- 
pative hydrodynamic calculations unless the proper time- 
dependence of the shear viscosity, which we have intro- 
duced in this paper, is taken into account. 
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